Particle-wave duality in quantum tunneling of a bright soliton 
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One of the most fundamental difference between classical and quantum mechanics is observed in 
the particle tunneling through a localized potential: the former predicts a discontinuous transmis- 
sion coefficient (T) as a function in incident velocity between one (complete penetration) and zero 
(complete reflection), while the later always changes smoothly as a wave nature. Here we report a 
systematic study of the quantum tunneling property for a bright soliton in ultracold atoms, which 
behaves as a classical particle (matter wave) in the limit of small (large) incident velocity. In the 
intermediate regime, the classical and quantum properties are combined via a finite (but not full) 
discontinuity in the tunneling transmission coefficient. We demonstrate that the formation of a 
localized bound state is essential to describe such inelastic collisions, showing a nontrivial nonlinear 
effect on the quantum transportation of a bright soliton. 

PACS numbers: 75.50.Pp, 75.30.Et, 72.25.Rb, 75.70.Cn 
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The interplay between the interaction effect and the 
disorder potential has long been an interesting subject 
in condensed matter physics, from Anderson localiza- 
tion in the noninteracting limit ll| to the Bose glass 
in the strongly interacting region [2[. Similar transport 
problem can also be investigated in the system of ul- 
tracold atoms, where Bose-Einstein Condensates (BECs) 
are demonstrated to unify concepts in classical and quan- 
tum physics at a macroscopic scale Fermionic or 
bosonic particles with a tunable interaction strength can 
be studied in a well-controlled quasi-disordered potential 
In this context, solitons, localized wavepackets un- 
dergoing confinement owing to nonlinear effects be- 
come an ideal representative for the investigation in a 
macroscopic scale of the wave-particle duality which is 
one of the fundamental pillars in modern physics Q • For 
example, a bright soliton (BS) resemble a classical par- 
ticle in their collision properties Q, and should have a 
complete penetration or reflection predicted by classical 
mechanics, see FigQJa); on the other hand, due to the 
underlying matter wave nature, a soliton should always 
reveal partial penetration and reflection as predicted in 
the quantum mechanics, see FigUJb). Therefore, it is 
of both interest and fundamental importance to study 
how the nonlincarity (i.e., the interaction effect between 
bosonic atoms) can modify the quantum transportation 
properties of BS and the transition between these two 
regimes [jj-flll|. 

Apart from the existing literature on optical spa- 
tial solitons fl2l - [l6j ]. in this paper we investigate quan- 
tum tunneling properties for a BS in one-dimensional 
(ID) bosonic atoms with different transport velocities. 
Through a systematic numerical simulation, we find that 
BS is like a classical particle(FigOTa)) when the incident 
velocity is small compared to the interaction energy and 



when the potential is repulsive, while it behaves as an or- 
dinary matter wave in the other limit and is independent 
of the sign of the local potential (FigQJb)). In the inter- 
mediate regime, the nature of particle- wave dualism from 
BS shows a discontinuity in the transmission coefficient 
(T) as a function of the incident velocity, while the ampli- 
tude of the discontinuity is less than one, as required by a 
true classical particle. We numerically calculate the full 
phase diagram in such a crossover regime, and observe 
a qualitative difference in the scattering process between 
a potential barrier and a potential well: the latter case 
is an inelastic scattering due to the appearance of a lo- 
calized bound state, sec FigJIJc). Semi-analytical curves 
for such a border are derived both for potential barriers 
as well as potential wells. The dual nature in quantum 
tunneling of BS elucidated in this work should be ready 
to be observed in the system of ultracold atoms as well 
as in the dielectric material with electromagnetic waves. 

Here, we consider the dynamics of a weakly inter- 
acting BEC at zero temperature, which can be well- 
approximated by the Gross-Pitaevskii equation, referred 
also as the nonlinear Schrodinger equation (NLSE) (l7l . 
GJ, 
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where the particle mass m and Ti are both set to 1, ^S(x, t) 
represents the condensate wavefunction, g measures the 
inter-particle interaction, and V(x) = Vq5{x) indicates 
a defect potential. When the interaction is attractive, 
g < 0, a stable bright soliton is supported in a uniform 
system with the solution 
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FIG. 1: (Color online) Schematic plots for different tunneling 
dynamics through a local potential, V(x). Vi, vt, and v r are 
the incident, transmission, and reflection velocities, respec- 
tively, (a) Classical picture with either a total transmission 
(T = 1) or total reflection (T = 0). (b) Quantum mechanical 
picture with a partial transmission (0 < T < 1). The no- 
tations, ^i^t, and $ r represent the incident, transmitted, 
and reflected wavefunctions. (c) Inelastic scattering process 
of a BS through a potential well (Vo < 0), where a localized 
bound state, appears after scattering. 



where the center of the wavepacket is denoted by x c , and 
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-Et, with the total energy E = vf /2+fx and 



the chemical potential fi = — /3 2 /2, respectively. The ve- 
locity for BS is characterized by Vi. By taking j3 = \g\/2 
to have a unit normalization, i.e., J_ CX} \^bs(x, t)\ 2 dx = 
2/3/ 1 g\ = 1, the only two independent parameters left are 
the normalized potential strength, Vb = Vo/\g\, and the 
normalized initial velocity, Vi = Vi/\g\, which defines our 
parameter space. In the following, we consider the trans- 
portation process when such BS wavepacket is generated 
at t = 0, centered at x c — > —00, and then propagates 
along the positive x-axis with an initial velocity v i . This 
wavepacket then scatters the local defect V(x) at the po- 
sition x = 0, resulting in possible transmitted, reflected, 
and localized wave functions after a certain time mea- 
sured. 

First of all, the calculated transmission coefficient (T) 
as a function of Vi and Vo > for a repulsive potential is 
illustrated in FigJ5fa) by directly solving Eq.(I) numer- 
ically. Here, T is defined as = J a °° lim^oo \^(x,t)\ 2 dx 
with a small value a > to exclude the contribution 
from any possible localized bound states. As one can see 
from Figj^a), in the region with a small value of Vi and 
I Vq I , there exists a line that characterizes the discontinu- 
ity in the transmission coefficient. The existence of such 
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FIG. 2: (Color online) Contour plot of the transmission coeffi- 
cient as a function of incident soliton velocity Vi and potential 
strength \Vo\. Cases of potential barrier Vo > and potential 
well Vb < are separately shown in (a) and (b). White dashed 
lines are plotted from Eq.([4]) and Eq.© for the corresponding 
analytical results (see the text). 



a discontinuity certainly reflects the particle nature of 
BS, i.e. totally transmitted (T = I) or totally reflected 
(T = 0) as shown in Fig. [IJa). However, this line of 
border for the particle nature breaks down at a critical 
point in the parameter space where the incident veloc- 
ity and corresponding potential strength are denoted by 
v* ~ 0.8 and \V * \ ~ 1.2. Beyond these values, instead of 
a disrupt change, the contours of transmission coefficient 
T changes continuously as a regular matter wave. 

On the contrary, in Fig(2Jb), the tunneling properties 
are different for a BS through an attractive potential 
(Vb < 0, a potential well). Qualitatively speaking, we 
have a similar "phase diagram" as the case of a potential 
barrier, but now the phase boundary to illustrate the par- 
ticle nature becomes a nearly linear one. As it would be 
demonstrated later, at a small value of Vq, such a univer- 
sal border, independent from any additional parameters, 
comes from the existence of a localized bound state. The 
formation of this localized bound state screens the poten- 
tial well and results in extra interactions and nonlincarity 
on the quantum tunneling of BS. 

In order to give a deeper understanding of these nu- 
merical results, we consider the limit of a weak nonlinear 
interaction for the first step. When the interaction effect 
is much smaller than the kinetic energy and potential 
energy, the corresponding transmission coefficient should 
be similar to that in the standard quantum mechanics 
textbook 0|, i.e., 
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As expected for a characteristic wave nature, the trans- 
mission coefficient T is always continuous and indepen- 
dent of the sign of potential strength, Vq. We note that 
above results are true for the incident wave as either a 
soliton wavepacket or a plane wave. In such a scenario, 
the incident BS can be easily distorted by the local po- 
tential due to that the nonlincarity is too weak to support 
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the original soliton solution, resulting in a lots of disper- 
sive radiations in the transmitted or reflected waves [19( . 

On the other hand, in the limit of a weak and repul- 
sive potential along with a small velocity, i.e., the strong 
interaction limit, we can safely assume that the propa- 
gating soliton is not affected by the potential. Hence, one 
can use the center position of BS, x c (t), to describe the 
whole transportation process if there is no bound state 
generated during the scattering process. In this limit, one 
can rigorously show that the dynamics of x c (t) behaves 
like a classical particle moving effectively in a conserva- 
tive potential, V c s{x), which is just a convolution of the 
local potential with the soliton wavefunction (20| . i.e., 

V e g{x) = ^jVo sech 2 (x). Therefore, the corresponding 
"conservation law" for the total energy is found to be, 



analytic form can be written as following: 
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where v(t) = -X dx J t is the defined particle velocity for 
BS. As a result, the border across a total reflection and 
a total transmission can be defined by taking v(t) = 
and x c = as the boundary condition, along with the 
initial condition Xi — > — oo. Then, we obtain the relation 
V /\g\ = (M//3) 2 (^/M) 2 = 4K/|<?|) 2 , which is depicted 
as the white dashed line in Figl^a). 

However, we know that above semi-classical approach 
used for Eq.Q fails when Vb is larger than a critical 
value, denoted as Vq* due to the breakdown of taking 
the BS as a classical particle. This critical value can be 
estimated as following: as the center of BS arrives the 
location of a potential, it gains a local potential energy, 
Vol^O)! 2 , which cannot be larger than the absolute value 
of the chemical potential, \n\ = f3 2 /2 = |g| 2 /8, in order to 
keep the soliton description valid. From Eq.@, we have 
the critical value for the breakdown as V^/lffl = 0-5. In 
FigEfa), the agreement between the analytical curve de- 
fined by Eq.((4|) and our direct numerical simulations is 
both qualitatively and quantitatively. When Vb is close 
to Vq, the BS is in the brink of collapsing, then the re- 
flection part as well as the non-soliton radiation become 
not negligible. Discrepancies from the above analytical 
results are therefore expected. 

Now we come to the potential well, which should have 
similar results as the potential barrier in the limit of large 
Vb and v% (the wave nature in the weak interaction limit). 
However, in the regime of small Vb and Vi, the wavepackct 
description used above fails for the lack in the consid- 
eration of possibility localized bound states supported 
in an attractive interaction. The appearance of a local- 
ized bound state indicates extra inelastic scattering and, 
therefore, changes the resulting tunneling amplitude dra- 
matically, as compared to the case of a potential barrier. 
The bound state wavefunction for a localized potential 
has been well-studied in the literature [21- 24|, and its 
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where (3b measures the slope and amplitude of the bound 
state, and Xb = tanh _1 (|Vb|//3i,) is the shift of effective 
peak position from the potential center. It is easy to see 
that the bound state wavefunction is composed of two 
soliton-likc solutions, but with different center positions 
and (3. By matching the discontinuity in the wavefunc- 
tion slopes with the potential strength, for a given renor- 
malization of the bound state, i.e., fib is fixed, such bound 
states exist only in a weak potential limit and disappears 
when | Vb| > fib- Such an anti-intuitive result originates 
from the fact that the maximum slope of BS is limited 
by its renormalization due to the nonlinear (interaction) 
effect. Although the localized bound state also exists in a 
repulsive potential defect, it cannot be easily produced in 
the tunneling process due to the mismatch in the bound- 
ary conditions. Therefore it docs not affect the tunneling 
property as we discussed above. 

Inspired by the numerical simulations of the tunneling 
process not shown her), we consider the following sim- 
plified picture of the tunneling dynamics in the presence 
of a bound state, i.e., inelastic scattering. To derive an 
analytical formula for the border when a soliton scatters 
by a potential well, we assume: (i) the reflected wave- 
function is negligible, (ii) the bound state appears after 
the scattering, and (iii) the transmitted wave also has a 
soliton profile. Since both the soliton solution and the 
localized bound state are governed by two parameters, 
P and v, as shown in Eq.([2]), the relevant parameters 
to describe a soliton tunneling are therefore: {(3i,Vi) for 
the incident soliton, {(3t,vt) for the transmitted one, and 
{/3b, Vb = 0) for the localized bound state. Since Vi is 
given and (3i = \g\/2 is required for the initial unit nor- 
malization, now we have only three parameters to be de- 
termined: /3b, f3 t , and v t . 

Instead of matching the boundary of wavefunctions 
during the scattering process, we use the conservation 
laws to extract the three unknown parameters above in a 
more general method. Based on the above three assump- 
tions, we can write down three equations: 
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where we have expressed the transmission coefficient 
2/3 t /\g\ as T via assumption (iii). Equations (JU) and 
represent the conservations of total probability and total 
energy by including both the localized bound state and 
the transmitted soliton. Equation © can be understood 
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as the conservation of current density due to the change 
of soliton amplitude. We note that |Vq| < A, in order to 
have T > in Eq.©. By eliminating the other two vari- 
ables, one can obtain T as a function of the normalized 
potential strength Vq = \ Vo/g\ and the normalized initial 
velocity Vi = Vi/\g\, 

_ Vo(V + 1) T Jvf(Vo + I) 2 - (4^o + 35?) 

T — , (9) 

2V + 3/2 

where the "+" solution is physically invalid. From 
Eq.(|9]), we find several interesting properties in the tun- 
neling of BS through a potential well. First of all, above 
solution for the transmission coefficient is real only when 
| V 1 2 ( | Vb I + 1 ) 2 - (4w 2 1 V 1 + 3u| ) > , or when v { is smaller 
than a critical velocity, w c (|Vb|), i.e., 

vc = \v /g\ (\Vo/g\ + 1) (w) 
\g\ ~ ^3 + A\V /g\ ' 

where Vq and Vi has to be bounded by requiring T < 
1. In FigfS^b), the border for the discontinuity in the 
transmission coefficient is compared with our analytical 
formula and direct numerical simulations, which results 
in very good agreement. More importantly, we find that 
the critical velocity defined in Eq. (fT0|) . i.e., the curve 
for the discontinuous transmission coefficient, shows a 
rather straight line (although is not exact), instead of a 
parabolic one for the potential barrier. The major reason 
is due to the existence of a localized bound state in this 
inelastic scattering process. 

Although our simplified theory does predict the loca- 
tions of transition, the transmission coefficient T we ob- 
tained in Eq.® is slightly less than the values obtained 
by direct numerical simulations. This is not surpris- 
ing because we have neglected the radiation parts (non- 
soliton waves) in the transmitted waves, which cannot be 
captured in such a simple theory. Last but not least, we 
restate that the estimations we made here for the tunnel- 
ing dynamics of a potential well cannot be applied to the 
potential barrier, because the bound state wavefunction 
supported by a potential barrier has to a double-humped 
one in the profile due to a mismatching boundary con- 
dition. On the side of wave nature, the transmission 
coefficient is identical both for potential well and poten- 
tial barrier when one considers a regular (noninteracting) 
matter wave tunneling, but on the side of particle nature, 
it turns out to be very different in the classical particle 
picture. This also reflects the non-trivial effect of inter- 
action (nonlinearity) in the soliton scattering problem. 

In conclusion, we have systematically studied the tun- 
neling of a bright soliton subjecting to a localized poten- 
tial defect. By performing direct numerical simulations, 
we obtain a full phase diagram of the transmission co- 
efficient in terms of the incident velocity and potential 
strength. Our results show an fundamentally important 



transport property, which can be easily observed in ultra- 
cold atoms, nonlinear optics, or even soft-matter systems. 
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